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ABSTRACT 



If a light gaugino sector exists, then the lightest supersymmetric particle (LSP) has a chance 

of being pair-produced in rare B-decays. As a consequence of neutral flavour violation in most 

><' 

supersymmetric models, such decays can occur at the tree-level and reinforce the channels 

B — > K(K*) + invisible. We discuss how a study of such decay spectra in B-factories can 

help us either identify or exclude a light LSP. 
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Since supersymmetry (SUSY) |T[ is one of the most appealing prospects in our quest for 
physics beyond the standard model, all possible ways to uncover (or rule out) its existence 
are worthy of attention, even if it be in relatively inprobable corners. Though the lower 
bounds on squark and gluino masses as announced by the CDF collaboration are 126 GeV 
and 141 GeV respectively 0, they are based on certain assumptions about superparticle 
decays and SUSY parameters. Relaxation of such limits cannot thus be altogether excluded 
||. In particular, there is a persistent, although controversial, claim that a window in the 
parameter space with a light gluino is still open @, |5| . If that indeed be the case , then a 
gluino in the mass range of, say, 2.5-5 GeV will cause a squark to decay directly into it. The 
gluino will subsequently decay into the lightest supersymmetric particle (LSP), supposedly 
even lighter in an R-parity conserving SUSY. As a result the missing transverse momenta 
(/Pt) associated with the LSP are considerably degraded and cannot survive the /Pt-cuIs 
employed in conventional SUSY searches. Therefore , this scenario also relaxes the squark 
mass limits M. Further motivations for a light gluino have come from the observation [B] that 
the values of the strong coupling a s at low and high energies are in better agreement with 
theory if a low-mass, electrically neutral coloured fermion is present. Attempts have been 
made to demonstrate the plausibility of such a scenario by proposing SUSY models where the 
gauginos acquire their masses radiatively f7|. On the phenomenological side, recent studies 
include the implications of a light, long-lived gluino in strong processes ||, constraints on 
light gluinos from electroweak precision tests |9j and those coming from radiative b-decay 

Here we focus on such a scenario from another angle; namely, we try to extract signatures 



of the light LSP from the kinematics of rare B-meson decays |ll|. We assume the LSP to be 
the lightest neutralino (x?) in this case. Then the flavour-changing neutral current (FCNC) 
decays B — > K(K*)x\x\ will give the same observable final states as B — > K{K*)uV. 
However, the decay spectra for the SUSY and standard model (SM) processes are going to 
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have different shapes. The net observed variation of dT /cIEk(k*) with the K(K*)-energy will 
be a result of superposition of the two types of final states, leading to a distribution with a 
kink. The position of the kink and the distortion to the spectrum relative to the purely SM 
case depends on the mass of the LSP. Some earlier studies had suggested the application of 
a similar effect in Kaon decay to unmask a massive tau- neutrino or a photino [|l!| . We have 
also recently discussed the implications of such spectral distortion in the context of decays 
of a heavy Higgs boson jn| . 



There are several advantages in looking for a light LSP in B-decays. First, as we shall 
see, an LSP in the mass range 0.5MeV — lGeV produces a kink in a conspicuous part of 
the spectrum, thereby making the distortion rather obvious. If we explore the window in 
rrig = 3 — 5GeV, then the LSP is likely to lie in the above range. Secondly, B-decays are 
dominated by short- distance physics, so that quark-level estimates are reliable. And finally, 
with a number of B-factories being designed for the near future, the prospect of producing 
10 7 ~ 8 BB-p&irs are realistic ||14| . Thus one may aspire to have a sufficient number of B's at 
the threshold so that their rare decays can be under scrutiny in the rest frame. 

Both the decays B — ► KxiXi an d B — ► K*XiXi are driven by the quark-level process 
b — ► sxiXi- in the supersymmetric standard model, there is the interesting posibility of 
tree-level flavour violation in quark-squark-neutralino (or quark-squark-gluino) interactions 



TS| . This can occur because the quark and squark mass matrices are not simultaneously 
diagonal. In a basis where the 3x3 down-quark mass matrix is diagonal, the 6x6 down- 
squark mass matrix is given by 
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where m s, m>u are the diagonal down-and up-quark mass matrix respectively, and K is the 



Kobayashi-Maskawa matrix, ttil and m^ are the flavour-blind SUSY breaking parameters 
that set the scale of squark masses. We can put tul = ttir without any loss of generality. 
The term proportional to m\ occurs due to radiative corrections induced by Yukawa coupling 
with charged Higgsinos. Such a term is particularly important in a model based on N=l 
supergravity (SUGRA), as the mass parameters evolve from the high SUGRA breaking scale 
to the scale of the residual global SUSY breaking at a lower energy. Because of this term, 
ra| cannot be simultaneously diagonal with m\ As a result, quark-squark-neutralino (or 
gluino) interactions can violate flavour. The flavour mixing in the down sector is controlled 
essentially by the top-quark mass, the Kobayashi-Maskawa matrix and the parameter cq. 
The calculation of cq is model-dependent; we shall treat it here as a phenomenologcal input 
that needs to be restricted by SUSY contributions to various FCNC processes. 

In the above framework (where we have also neglected the mixing between left-and right- 
squarks), the quark-squark-neutralino (q — q — Xi,i = 1—4) coupling in the down sector 



is 
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q k + h.c. (2) 



where Tj k is the (jk)-th element of the unitary matrix that diagonalises the upper 3x3 block 
of m| in equation(l). N is the neutralino mixing matrix, and T%j the third component of the 
the isospin of the j-th flavour. In the absence of left-right mixing, only the first term in (2) 
is relevant for flavour-violating interactions. 

For our calculations, the element T 2 3 will be important. We have used m t = 170GeV. 
For such a top-quark mass, the third term in the upper-left block of m| is important from 
the viewpoint of diagonalisation, so that for a not-too-small value of Co, the elements of T are 



close to those of K in magnitude. We have parametized T 2 3 by writing r 23 = cK 2 -a- Various 
values of Cq and and the corresponding values of c are given in table 1. 

The two tree-level graphs (plus those with the momenta of the neutralinos interchanged) 
shown in figure 1 contribute to b — ► sxiXv F° r numerical calculations, we confine ourselves 
to a situation where the gluino mass is in the range 3-5 GeV. In such a case, the other 



parameters in the SUSY sector have to be compatible with the LEP-I results |16] . Under such 
circumstances, it is straightforward to verify that the LSP is dominated overwhelmingly by 
the photino. This is easily demonstrated if, for example, we adhere to a scenario inspired by 
Grand Unified Theories (GUT) |17|| . In such cases the only independent inputs apart from the 
gluino mass are u, the Higgsino mass parameter, and tan/3, the ratio of the two scalar vacuum 
expectation values. In the above situation one gets confined to 50GeV < /i < lOOGeV and 
1.0 < tan/3 < 1.5. Diagonalisation of the neutralino mass matrix with such inputs reveals the 
complete dominance of the photino state in the LSP. Therefore, the interaction (2) can be 
used in the X? = 7 limit for our purpose [1(b)]. The amplitude for b(p ) — > s(p 3 )xi(p2)Xi(Pi) 
can be expressed as 
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where nig is the assumed common mass of the s and b squarks. In order to factor out the 
hadronic part of the amplitude, one needs to do a Fierz transformation on (3), which gives 
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77 72 ^ s (Ps) 7 M (1 - 75) KPo) u (Pi) 7 m (1 + 75) v(p 2 ) - 

[{Po ~ Pi) ~ mf\ 

-2 s(p 3 ) Y (1 - 7s) 6(po) «(P2) 7m (1 + 75) v{p x ) (4) 
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To proceed further, one has to obtain the hadronic matrix elements for the above quark 
current. Using a common parametrization for the matrix elements for rare B-decays ||18| , 
one can write 

(K(p 3 )\s^b\B(p )) = f + (q 2 )(p 3 +p ) fi + f4q 2 )q, , (5) 

(K(p 3 )\s 1 , l5 b\B(p ))= (6) 

(K*(p 3 )\ S J, (1 - 75 ) b\ B(p )} = i e, uaP e"(p 3 ) (p + p 3 ) a q^V(q 2 ) - 
Cm(Ps) [ m l - m2 K* } Mq 2 ) -(eg) (Po+Ps)^ A + (q 2 ) - (e . q) g M A_(q 2 ) (7) 

where g M = (p — p 3 ) and e /i (p 3 ) is the polarization vector for the K*. Our results are based 
upon numerical values of the various form-factors (and pole fits for their momentum-transfer 
dependence) obtained from the relativistic quark model of reference [|n| . These form-factors 
have been computed in the literature using other models, too |2U[|; the question of matching 



such calculations with the general relationship predicted by heavy quark effective theories 
have also been discussed [OTJ. However, the uncertainties due to model-dependence do not 



affcet the general features of our results. It should also be noted that we have not taken 
QCD corrections into account here. Though such corrections moderately alter the decay 
rates [22|, the key featurs of the spectral distortions should not be affected, since at the 



lowest order electroweak level, the SUSY and standard model effective interactions have the 
same operator structure. 

The differential decay rates of our interest are given by 

dT 1 r Ei(max) r 



(IEk (K*) 647T 3 VTiB JEi(min) 

where \M\ is the squared matrix element and 
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Ei (max) 



(m B - E K (#*) + p 3 \/l - 4m£o/<? 
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Ei(min) = - 



[m B - E K ( K *) - p 3 Jl - 4m 2 x o/q 
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q 2 = m\ (Jf *j + m 2 B - 2 m B E K (K *) (11a) 

P3= E K(K*)- m K(K*) ( llb ) 

^k - (k*) is the i^ (if*) energy in the rest frame of decaying B and its kinematically allowed 
range is 

m K (k*) < E K (k*) < (m% + m 2 K {Kir) - 4m x o)/2m B (lie) 

Further, one has to add the rates for the SUSY process with that for HB — > K(K*)viVi 



which occurs via triangle as well as box diagrams |23 . 

The numerical results are shown in figures 2-5. For the K* final states, the polarizations 
have been summed over. The fact that the standard model final states consist of vv pairs 
explains why the distributions are not vitiated by peaks due to the J/ip and tpi resonances. 
We have drawn the graphs for rriq = lOOGeV^ which is easily allowed in this scenario. Two 
sets of graphs each for the K and K* final states are presented in order to demostrate the 
dependence of the rates on the parameter cq. Evidently, with even quite conservative choices 
for Co one can notice distortions to the spectrum over a considerable region of the parameter 
space. The effect becomes less and less obvious with increasing squark mass, and is barely 
perceptible for B — > KxiXi with cq ~ .01, rriq = 500GeV. Also, the response to a variation 
in the mass of the LSP in the region 0.4 — lGeV is manifest. A few hundred events in a 
B-factory should suffice to explore this kind of a distortion. 

It is to be noted that while the differential decay rate for E(S — ► Ki>iVl) increasess 
monotonically with Ek, it dips after an initial rise in the case of E(S — > K*viUj). This is 



because the transverse component of K* has an important role in the rate. For it, the upper 
end of the phase space (corresponding to the maximum allowed value of Ex* ) corresponds to 
a configuration that is disfavoured by helicity conservation. The corresponding distortions 
to the spectrum caused by the LSP are less conspicuous than in the case of decays into a K, 
although the overall rate for the former is higher by about an order of magnitude. 

The graphs also clearly portray the fact that with light LSP's the total rates for both 
the decays can be jacked up by as much as an order of magnitude due to presence of a light 
LSP. At least over a certain amount of parameter space such an overall rate enhancement 
can be taken as a positive signal of such an LSP. 

To conclude, if 10 7-8 SB-pairs are produced in a B-factory per year, then the tree- 
level flavour violating interaction in SUSY models strongly affect the decay patterns in 
B — > K + nothing and B — >■ K* + nothing , provided that a light LSP is present. This can 
be, in a somewhat model-dependent way, a pointer to the light gluino as well. Therefore, 
a detailed investigation of the above types of decays in B-factory experiments are going to 
help one in constraining the light sparticle scenario to a large extent. 
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Table caption 

Table 1: Different values of the parameter Cq and the corresponding values of c, for rriq 
100 GeV. 
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Figure Captions 



Figure 1: 
The tree-level contributions to b — ► sXiXi- in addition there will be crossed diagrams where 
the four-momenta of the LSP's are interchanged. 

Figure 2: 
The differential decay rates for B — >K + nothing for niq = lOOGeV, c = 0.1. The solid, 
dotted and short-dashed curves correspond to three LSP masses expressed in GeV. The long- 
dashed curve below is for the purely standard model case with three massless neutrinos. 

Figure 3: 

Same as figure 2, but with c = 0.5. 

Figure 4: 
The differential decay rates for B — >K* + nothing, with the same choice of parameters as 
in figure 2. 

Figure 5: 
Same as in figure 4, but with c=0.5. 
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